UPPER BOUNDS FOR THE NUMBER OF LIMIT CYCLES OF 
SOME PLANAR POLYNOMIAL DIFFERENTIAL SYSTEMS 



ARMENGOL GASULL AND HECTOR GIACOMINI 



Abstract. We give an effective method for controlling the maximum number 
of limit cycles of some planar polynomial systems. It is based on a suitable choice 
of a Dulac function and the application of the well-known Bendixson-Dulac 
Criterion for multiple connected regions. The key point is a new approach to 
control the sign of the functions involved in the criterion. The method is applied 
to several examples. 



1. Main result 

One of the few general methods that allows to give upper bounds for the number 
of limit cycles of planar differential systems 

x = P(x,y), y = Q(x,y) 

is the use of Dulac functions in multiple connected regions, see [21 El ® El [HI E] ■ 
Recall that the primer idea is that when the function div(P, Q) does not vanish 
on a simply connected region U C M 2 , then the above differential system has 
no periodic orbit totally contained in U. We state the general Bendixson-Dulac 
Criterion in next section, see Theorem 12. II The main difficulty for practical uses of 
this result is that it is needed to find a (Dulac) function g such that div(gP, gQ) 
does not vanish on a suitable set. This paper gives a quite general result for 
polynomial differential systems, see Theorem El Its proof is based on a "good" 
choice of a Dulac function. As we will see in Section [3] this result provides a 
constructive way for giving upper and lower bounds for the number of limit cycles 
of a large class of planar polynomial systems. 

Given a polynomial p(s) G M[s], we will say that the couple (k,w(r)) e M + x]R[r] 
is a Dulac pair of p(s) if 

Pk,w(r) := rp(r 2 )w' (r) — 2k (p(r 2 ) + r 2 p'(r 2 )) w(r) < for all r > 0. 

As we will see in Lemma ET] and in the proof of Proposition [221 the above inequal- 
ity implies that the function |w(r)| _1 / fe is a Dulac function, in any of the connected 
components of M 2 \ {w(r) = 0}, for the systems that in polar coordinates, r and 
8, writes as f = rp(r 2 ), 8 = p(r 2 ). Here p is any arbitrary polynomial. We remark 
that there is no need to introduce this Dulac pair to study the phase portrait of 
this system. For instance their limit cycles are given by the positive zeros of p 
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that are not zeros of p. Nevertheless, as we will see in next theorem, Dulac pairs 
are useful to study more general systems. 

Theorem A. Consider the polynomial differential system 

x = P(x,y), y = Q(x,y), (1.1) 

where P and Q are real polynomials of degree n, satisfying that P(0, 0) = Q(0, 0) = 
0. In polar coordinates it writes as 

r = R(r, 9) := P(r cos 9, r sin 9) cos 9 + Q(r cos 9, r sin 9) sin 9, 
9 = Q(r,9):=- (q{t cos 9, r sin 9) cos 9 — P{r cos 9, r sin 9) sin 9^j . 
Define the polynomial 

1 f 2n 
p( r 2 ) ■= / R(r,9)d9 

and assume that p(s) has the Dulac pair (k,w(r)), where w is a polynomial of 
degree d. For these k and w consider the function 

\ n/ „, „ x , fdR{r,9) dQ{r,9) R(r,9)\ . . 
Mr, 9, k,w): = R(r, 9)w'(r) - k — ^-i- + — ^-^ + v ; w(r) 

\ or 69 r J 

n+d-l 

=: M t {9,k,w)r\ (1.2) 

i=l 

and, for any i > 1, Zet rrii(k,w) be such that max^o^d Mj(0, /c, w) < rrii(k,w). 
Let m + be the number of non-negative roots ofw. Then, if the polynomial 

n+d-l 
i=l 

is negative for all r e (0, oo), system (II .ip has at most m + limit cycles and all of 
them are hyperbolic. 

Next remarks collect some comments on the above result. 

Remark 1.1. (i) The above theorem gives a constructive way of testing when a 
Dulac function of a system of the form r = rp(r 2 ), 9 = p(r 2 ) is also suitable for 
system (11. lj) . 

(ii) In Proposition \2.5\ we will prove that there are many polynomials p for 
which a Dulac pair exists. Moreover, in Section^ we will show many systems of 
the form (II. 1ft for which the above theorem can be applied. 

(Hi) Remark \2. 81 proves that the upper bound given in the theorem can not be 
improved. 

Remark 1.2. (i) It is not difficult to see that any differential system fulfilling the 
hypotheses of Theorem\M has at least m + — 2 — C limit cycles, where C > can be 
computed as follows: Let < n < • • • < r m + be the ordered non-negative roots of 
w(r). Define the rings Ri = {(x, y) : r« < ^ x 2 + y 2 < r i+ i} } i = 1, 2, . . . , m + — 1. 
Then C is the number of rings among R±, R2, . . . , R m +-i containing critical points 
of system (II. ip . 
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(ii) From the proof of the theorem it follows that there is at most one limit 
cycle in each of the sets R\, R 2 , ■ ■ ■ , R m +-i, R m +, where R m + is the unbounded ring 
R m + = {(x,y) : r m + < \J x 2 + y 2 } and that (when r\ > 0) there is no limit cycle 
in the disc {(x, y) : \J x 2 + y 2 < r±}. Moreover in the rings R 2 , R3, ■ ■ ■ , R m +-i not 
containing critical points of the system there exists always a limit cycle and its 
stability is given by the sign of w on them. The limit cycle on R\ only can exist 
when ri > or when T\ — and the sign of w on this set does not coincide with 
the stability of the origin. 

Let us see how the theorem works in a concrete example. We will prove that 
the system 

49 26 1 

x = — y + Ax x xy H — xy + x + 2x y + xy = P(x,y), 

10 5 5 

23 12 
y = x + Ay- —x 2 y - 5y 3 - -xy 3 - —y A + x 4 y + 2x 2 y 3 + y 5 = Q{x, y) (1.3) 
5 5 15 

has exactly two limit cycles. Examples with parameters will be studied in Sec- 
tional Following the notation of the theorem we get that p(s) = A — 79s/16 + s 2 . 
Moreover, by using Proposition 12.51 we take w(r) = r 2 p'(r 2 ) = r 2 (2r 2 — 79/16). 
Finally we choose k = 4/5. We have 

79 2 1287 4 237 , 8 s 

Ph m [r) = r r H r r . 

^ k ' wy 1 10 128 40 5 

By using Sturm's rule it is easy to prove that Pk, w (j) < for all r 7^ and so 
(k,w(r)) is a Dulac pair for p. Straightforward computations give that 

M(r, 6, k, w) = r 2 + f -25740 + cos(2fl) + 316 cos(40) J r 4 



10 2560 

79 

12 cos(#) + 3 cos(30) - 15 cos(50) + 46 sin(^) - 7 sin(3#) - 5 sin(50)) r 
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4800 



1 / 1 oo \ 

+ — U74 - — cos(20) - 8 cos(4#) J r 6 

+ (6 cos(^) + 9 cos(3#) - 15 cos(50) - 2 sin(0) + 9 sin(30) - 5 sin(50)) r 7 - ^r 8 . 



By using rough bounds like 

-46 < 6 cos(#) + 9 cos(3fl) - 15 cos(50) - 2 sin(0) + 9 sin(30) - 5 sin(50) < 46 
we get that 

B *-/ /i 7 \ t / \ 79 2 3459 4 869 5 1269 6 46 7 8 8 

M(r, 9, k, w) < ^kw(r) = r 2 r 4 H r 5 H r b H r 7 r 8 

v ' ' ' ; - ' v ; 10 400 600 200 75 5 

and it can be proved by using Sturm's rule that it is negative for all r > 0. 

Therefore, since w has two non-negative roots, m + = 2 and the system has at 

most two (hyperbolic) limit cycles. 

It is clear that the origin is an unstable focus. By taking the resultant of P 
and Q, for instance with respect to y, and applying Sturm's rule to this new 
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polynomial we prove that the origin is the unique critical point of the system. By 
studying the flow of the system on the circles {x 2 + y 2 = R 2 } for R big enough, 
and on {x 2 + y 2 = r 2 }, where r is the positive root of w, and using the above 
information we deduce that the two limit cycles actually exist. Indeed one of them 
is in P = {x 2 + y 2 < r 2 .} and is stable and the other surrounds T> and is unstable. 

In fact, we have presented the study of system (11. 3ft in this introduction because 
it shows at the same time the power and the limitations of the method. By one 
hand, for a concrete system which has no a priori structure, we can give the exact 
number of limit cycles. On the other hand, in all the cases that we have studied, for 
polynomials systems of degree 2n or In + 1 it only has been applicable for systems 
having at most n limit cycles and when all these limit cycles are hyperbolic and 
nested. 

Due to the extreme difficulty of the second part of the Hilbert's sixteenth prob- 
lem it is never easy to give explicit and realistic upper bounds for the number of 
limit cycles of a planar polynomial system. Moreover most results in the literature 
can only be applied to control the limit cycles of particular types of differential 
systems, namely Lienard systems, quadratic systems, cubic systems, systems with 
homogeneous nonlinearities, etc, giving usually upper bounds of 0, 1 or 2 limit cy- 
cles, see for instance [TJ EU [H] . The criterion proved in Theorem lAl is not subject 
to these restrictions. 

The paper is organized as follows: Section 2 contains some preliminary results 
and the proof of Theorem |A] In Section [3], we will apply this result to different 
families of differential systems providing in all of them explicit upper bounds for 
the number of limit cycles. In fact, we present two examples with limit cycles 
surrounding a unique critical point (see Examples 1 and 2), two situations with a 
limit cycle surrounding several critical points (see Examples 3 and 4) and finally a 
family of examples of differential systems showing that, in the set of systems of the 
form (11.11) . there are many open subsets of systems satisfying all the hypotheses 
of the theorem. 



2. Preliminary results and proof of Theorem A 
Let U C M 2 be an open set with smooth boundary and such that its fundamental 

£ 

group, 7Ti(W), is Z* • • • *Z or in other words having £ gaps. For short we will say 
that U is an £-punctured open set and we will denote by £(U) its number of gaps. 
Notice that with this notation, simply connected sets are 0-punctured sets and 
multiple connected sets are ^-punctured sets with £ > 1. 

The following result is a well known extension of the Bendixson-Dulac Criterion 
to ^-punctured sets. For a proof, see any of the papers [HI El E] • 

Theorem 2.1 (Extended Bendixson-Dulac Criterion). Consider the C 1 -differential 
system 



x = P(x,y), y = Q(x,y) 



(2.1) 
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and set X = (P,Q). Let U be an open £-punctured subset of M 2 with smooth 
boundary. Let g: U — > R be a C 1 -function such that 

M := divfoX) = d JL P+ d J- Q + ( ?L + ?Q )g = (V9 , x) +gdiv(x) 

does not change sign in U and vanishes only on a null measure Lebesgue set and 
such that {M = 0} fl {g = 0} does not contain periodic orbits of ^2.1\) . Then the 
maximum number of periodic orbits of Ii2.1\) contained inlA is £. Furthermore each 
one is a hyperbolic limit cycle that does not cut {g = 0} and its stability is given 
by the sign of gM over it. 

We will prove and use the following corollary of the above theorem. 

Proposition 2.2. Consider the C 1 - differential system 

x = P{x,y), y = Q{x,y) (2.2) 

and set X = (P,Q). Assume that there exist a positive real number k and a 
polynomial f(x,y) such that 

Mk := ^P+^-Q~ k(^ + = (V/, X) - kf div(X) (2.3) 

ox ay ox ay 

does not vanish in an open region with regular boundary Wet 2 . Let U\, U2,- ■ ■ , 
U m , be the connected components of W \ {/ = 0}. Then 

(i) The periodic orbits of system ( 12.21) contained in W never cut the set {/ = 0}. 

(ii) The maximum number of limit cycles of system (12.21) contained in each Uj, 
j = 1,2, ... ,m, is £(Uj), and all of them are hyperbolic. Moreover their 
stability is given by the sign of —fM^ on each region. 

(Hi) The maximum number of limit cycles of system (12. 2p in W is £(Ui) +£(^2) + 
• • • + £(U m ) and all them are hyperbolic. 

Proof. Since does not vanish on W we know that (V/, X)\ rj =0 } n >v does not 
vanish. Therefore the periodic orbits of (12.21) which are totally contained in W 
never cut the set {/ = 0} and (i) follows. To prove items (ii) and (iii) it suffices to 
apply Theorem 12. II to each one of the connected components W \ {/ = 0}. Note 
that, on each of them, the function g = |/| _1//fc is smooth and moreover 

div(^x) = zE^Z)|/|-i/*-i [(V /,x) - fc/div(x)] = zE^|/|-i/*-i Mfc . 

Therefore the upper bound for the number of limit cycles follows. 

Finally we prove the hyperbolicity of all the limit cycles. Let 7 = {(x(t), y(t)) : 
t G [0, T]} be one of them. We must prove that f Q div X(x(t),y(t)) dt ^ 0, see for 
instance Thm. 1.23]. Since 7 does not intersect the set {/ = 0} we have that 
over 7 

,. „ M t (V/.A') 

Hence 
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and its stability is given by the sign of —fM^ on the region where 7 lies. Therefore 
the result follows. □ 

To apply the above result, it will be useful to get the expression of the function 
Mfc given in the Proposition 12.21 in terms of the components of the differential 
system (12.21) . written in polar coordinates. 

Lemma 2.3. Let r = R(r,8), 6 = 0(r, #) be the expression of system (12. 2D in 
polar coordinates. Then the function (12. 3D writes as 

ox oy ox oy 

- ^L R + ^Iq - k (— + — + -) f 

dr 86 \dr 86 r J 

Lemma 2.4. Let p(s) be a real polynomial having all its roots real and simple. 
(1) For alls G R 

p(s)p"(s) - (p\s)) 2 < 0. 
(ii) For each fcel define the new polynomial 

q k {r) := 2r 4 (p{r 2 )p" (r 2 ) - k(p'(r 2 )) 2 ) + 2(1 - k)r 2 p{r 2 )p (r 2 ). 

Then there exists k, with \k — 1| small enough, such that qk{r) < for all r ^0. 

Proof. It is not restrictive to write p(s) = Y[ l j=i( s ~~ s j) "with S{ 7^ Sj for i 7^ j. 
From the equalities 



p(s)p"(s) - (p'(s)f 




item (i) follows. 

To prove (ii) notice that by item (i) 

gi(r) = 2r 4 (p(r 2 )p"(r 2 ) - (p'(r 2 )) 2 ) < for all 

Moreover qk(r) = qi(r) + (1 — k)q(r) for some polynomial gof degree 41 and 



q h (r) = 2coCi(l - k)r 2 + 2 (c 2 + 4c c 2 - 2k(c\ + c c 2 )) 



1 



+ ■■■ + 2{{l- k)l 2 -kl) r il , 

where p(s) = c + cis + c 2 s 2 H hs ; . Note that p(0)p" (0) - (p'(0)) 2 = 2c c 2 -cf < 

and thus, for \k — 1| small enough, (1 — k)l 2 — kl and c 2 + 4cqC 2 — 2k{c\ + C0C2) are 
both negative. Hence when cqCi = 0, taking \k — 1| small enough, we can always 
ensure that the sign of qkij) is negative in E\ {0}. When cqC\ 7^ 0, in addition we 
have to take k such that cqCi(1 — k) < 0. □ 
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Proposition 2.5. Let p(s) be a real polynomial. Then 

(i) Ifp(s) has all its roots real and simple by taking w(r) = r 2 p'(r 2 ) and some k 
near 1 then (k,w(r)) is a Dulac pair ofp(s). 

(ii) Ifp(s) has some real multiple positive root then it has no Dulac pair. 

(Hi) There are many polynomials p(s) with complex roots also having Dulac pairs. 

Proof, (i) Note that 

Pk,w(f) = rp{r 2 ) (r 2 j/(r 2 ))' — 2k (p(r 2 ) + r 2 p'(r 2 )) r 2 p'{r 2 ) = qk(r), 

where q^ is given in Lemma 12.41 By using this lemma the result follows. 

(ii) Let s* > be one of these roots. It is easy to see that for any couple k and 
w, the polynomial Pk, w ( r ) also vanishes at y/s*. Hence the result follows. 

(iii) We construct a class of polynomials for which the result holds. Let p(s) 
be a polynomial such that p(s) = Pi{s)p2{s) where pi(r 2 ) + r 2 Pi(r 2 ) > for all 
r G R and P2{s) = Yii=i( s ~ r f)i being rf different positive numbers. Consider 
w(r) = P2(f 2 ) and k — 1. Then 

p 1)W (r) = -2pj(r 2 ) {pi{r 2 ) + rVi(r 2 )) < for all r el 

Take now w e (r) = P2{r 2 ) + e 2 p' 2 (r 2 ). Then 

Pi,w £ ( r ) =PiA r ) +£ 2 W(r) 

for some polynomial W of the same degree that p\ >w . Moreover 

W(r l ) = -2r 2 p 1 (r 2 )(p' 2 (r 2 )) 2 <0 

for all i — 1, . . . ,j. Hence for |e| small enough the polynomial pi Ws (r) is negative 
for all r G R, as we wanted to prove. □ 

Remark 2.6. Notice that, given a polynomial p under suitable hypotheses, item 
(i) of Proposition HT51 provides a constructive way of finding Dulac pairs. On the 
other hand item (iii) provides a theoretical way to see that the same situation holds 
for other polynomials. Nevertheless, for a given polynomial p, even with complex 
roots, it is not difficult to find a Dulac pair. We give a couple of examples and 
some intuition of how we get them. 

(a) Forp(s) = (s - 2)(s - 4)(s 2 + 4)(s + 3) takew(r) = (r 2 -l)(r 2 -3) and k any 
of the values 1/2, 1 or 2. 

(b) For p(s) = —35 — 36s + 49s 2 /2 — 14s 3 /3 + s 4 /4, which has a unique positive 
root sq ~ 11.12 take k = 1 and w(r) = r 2 — a for any a for instance in the 
interval [5.5, 11]. 

By Remark li.il the existence of k > and w such that pk, w < for all r > 
implies that we can apply Proposition \2.S\ with f = w(r). Therefore each one of 
the hyperbolic limit cycles of system f = rp(r 2 ), 9 = p(r 2 ), which are given by the 
positive simple zeros of p, has to be contained in one of the 1-punctured regions 
o/R 2 \ {w = 0}. For instance in case (b), there is only one limit cycle r 2 ~ sq 
and we can try with a function w(r) = r 2 — a with a smaller than this value. We 
construct the function w of item (a) in a similar way. 
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Next lemma will be useful to get systems for which the hypotheses of Theorem 1X1 
hold and to prove that the upper bound given by the theorem is optimal. 

Lemma 2.7. Consider system 

x— xu{x 2 + y 2 ) — yv (x 2 + y 2 ), 

y— xv(x 2 + y 2 ) + yu[x 2 + y 2 ), (2.4) 

where u and v are arbitrary real polynomials. If u has all their roots real and 
simple then, taking w(r) = r 2 u'(r 2 ), the function $i jM (r) introduced in TheoremlA\ 
is negative for all r e (0, oo). 

Proof. Writing the system in polar coordinates we obtain 

f— R(r,9) — ru(r 2 ), 

0= G(r,6)= v{r 2 ). 

Clearly, from the above expression we get p = u, where p is the polynomial intro- 
duced in Theorem lAl Thus, taking k = 1 and w(r) = r 2 u'(r 2 ), we obtain 

M(r, 6>, 1, w) = R(r, 9) (r u (r )) - I — — 1 — 1 — J r u (r ) 

= ru ( r 2 ^ ^r 2 u'(r 2 ))' — ^(ru(r 2 ))' + u{r 2 ) \ r 2 u{r 2 ) 
= 2r 4 (u(r 2 )u"(r 2 ) - {u\r 2 )f^ = $^(r). 
By using Lemma [2.41 (i) the result follows. □ 

Remark 2.8. Notice that the above lemma implies that the upper bound given in 
Theorem [31 can not be improved. Consider in system (12 .4p a polynomial u such 
that all their roots are real and has no common roots with v. Then system (12. 4ft 
has as many limit cycles (indeed invariant circles) as number of positive roots, 
say m* . It is easy to take a polynomial u under the above hypotheses such that 
v! has exactly m* — 1 positive roots. Hence the number of non-negative roots of 
w(r) = r 2 p'(r 2 ) = r 2 u'(r 2 ) is m + = m*. By applying Theorem L"4l we get an upper 
bound of m + = m* limit cycles, which is indeed the actual number of limit cycles 
of the system. 

Proof of Theorem \M We apply Proposition l2.2l to system (II. ip with f(r, 8) = w(r) 
and the value k given in the statement of the Theorem. Then, by using Lemma 
we get that the expression of given in Proposition 12.21 is 

M k = R(r, Wr) - k + + ^) »(r) = M(r, 6, k, w ) 

where M(r, k, 8, w) is the function given in Theorem El 
Moreover, by hypothesis, we have 

n+d—l n+d—1 

M k = M{r,8,k,w) = Mi{9,k,wy < ^ m^k, w)r* = $ fc)W (r) < 

i=l i=l 
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for all r e (0, oo). Therefore, by Proposition ^. 2\ the maximum number of limit cy- 
cles can be bounded above by studying the topology of the connected components 
of the set W := M 2 \ {w(r) = 0}. Clearly it has m + connected components, all of 
them indeed 1-punctured when u>(0) = and it has m + + 1 connected components 
when w(0) ^ 0. Notice that in this later case one of them is simply connected 
and the other m + are 1-punctured sets. In any case, again by Proposition 12.21 
there is no limit cycle in the simply connected component and there is at most 
one limit cycle in each of the 1-punctured components, which is hyperbolic when 
it exists. Moreover, since < 0, its stability is given by the sign of w on each 
component. As it is already said in Remark 11.21 (h). it can be proved that the 
bounded 1-punctured components of W not having critical points of system (II. ip 
contain effectively a limit cycle. This result holds because each one of these rings 
is either positively or negatively invariant by the flow of the system. □ 

3. Examples 
3.1. Example 1. Consider the system 

x = x(l — (x 2 + y 2 )) — y(l + 2{x 2 + y 2 )) + axy + bxy 2 , 

y = x(l + 2{x 2 + y 2 )) + y(l - (x 2 + y 2 )) + cy 2 + dx 3 , (3.1) 

which has at the origin an unstable focus. Let us see that when b < 8 and a, c and 
d are such that 

* aACi(i (r) := -12 + (2\2a - c\ + 10|c - a\) r + (2b - 16 + 12|6| + 15|d|) r 2 < 

for all r > 0, then it has at most one limit cycle, which when exists is hyperbolic 
and stable. 

We apply Theorem IA1 Then p(s) = ((b — 8)s + 8)/8. By using Proposition 12.51 
we take w(r) = r 2 p'(r 2 ) = (b — 8)r 2 /8 and we choose k = 7/10. Then 

M(r, 6, k, w) = 3(& ~ 8) r 2 + — ( sin(0) + (c - a) S m(36)) r 3 
40 16 \ 5 / 

+ h -j^- ( 16 ~ 2b + b cos(2#) - cos(A9^j +d(2 sm(29) + sin(40))^ r 4 . 

Analogously that in the example given in the introduction, we have 

g ^ 

M(r,9,k,w) < —77-r 2 ty aM4 (r) < 
160 

for all r > and we are under the hypotheses of the theorem. Since m + = 1 the 
uniqueness of the limit cycle follows. It is easy to see that for a, b, c and d small 
enough the condition on ^ a ,b,c,d holds and the limit cycle exists. 

Note that for this system we can give a simple and explicit condition on the 
parameters of the system under which we can prove that there is at most one limit 
cycle. For instance the condition holds for a = c = 1 and —2b = 2d = 1. For these 
values we have also checked numerically that the limit cycle actually exists. 
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3.2. Example 2. Consider the system 

x = x(l — (x 2 + y 2 )){2 — (x 2 + y 2 )) — y + ax 2 y + bx 2 y 2 , 

y = x + y(l-(x 2 + y 2 )){2 - {x 2 + y 2 )) + cxy 2 . (3.2) 

If a, b and c are such that 

*a,&,c(r) := -10 + ~(|a| + |c|) + ~\b\r + (12 + \a\ + \c\) r 2 + |6|r 3 - 4r 4 < 

for all r > 0, then system (13.21) has at most two (hyperbolic) limit cycles. Moreover, 
when they exist, one is included in the disc V := {x 2 + y 2 < 3/2} and is stable 
and the other one is outside the disc and is unstable. 

The proof follows again by using Theorem S We take p(s) = 2 — 3s + s 2 , and 
by Proposition 12. 5\ we consider w(r) = r 2 p'(r 2 ) = r 2 (— 3 + 2r 2 ) and k — 1. Then 

M(r, 9, k, w) =- (-40 + a (6 sin(2^) - 3 sm(49)) + c (6 sin(2#) + 3 sin(40))) r 4 

+ h (2 cos(9) - 3 cos(3fl) + cos(5^)) r 5 
8 

b 

+ (12 + a sin(4#) - c sin(40)) r 6 - - (- cos(30) + cos(50)) r 7 - 4r 8 . 

Hence, for the values of the parameters considered, we have 

M{r,9,k,w) <r^ aAc {r) <0 

for all r > 0. Thus we can apply Theorem |X] with m + = 2, proving the existence 
of at most two (hyperbolic) limit cycles. 

For instance the condition on ^ a ,b,c holds for a = 1/8, b = 1/15 and c = 1/20. 
Moreover for these parameters it is not difficult to prove, by using resultants and 
the Sturm's rule, that the origin is the unique critical point, which is unstable. 
Finally, by studying the flow on {x 2 +y 2 = R 2 }, for R big enough, and on {x 2 +y 2 = 
3/2}, we prove the existence of both limit cycles. 

3.3. Example 3. Consider the system 

x = x(l- (x 2 + y 2 ))(2 - (x 2 + y 2 )) - y(l - (x 2 + y 2 )) + ax 4 + bx 2 y 2 , 

y = x(l- (x 2 + y 2 )) + y(l - (x 2 + y 2 )){2 - (x 2 + y 2 )). (3.3) 

If a and b are such that 

m ah (r) : = _10 + 27 \ a \+ 9 \ b \ r + 12r 2 + (2\a\ + |6|) r 3 - 4r 4 < 

for all r > 0, Theorem lAl will allow us to show that system ( 13. 3ft has at most two 
limit cycles. Moreover we will see that when there is no critical point outside the 
disc V := {x 2 + y 2 < 3/2} an unstable hyperbolic limit cycle always exists outside 
the disc V and it surrounds several critical points. 
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We take p, w and k as in the previous example. Then 
M(r,9, k, w) = -10r 4 + 12r 6 - 4r 8 

+ - (a (14 cos(#) + 3 cos(3#) - cos(5#)) + b (2 cos(fl) - 3 cos(3#) + cos(50))) r 5 
8 

+ - (a (-2 cos(#) + cos(3#) + cos(5#)) + b (cos(3#) - cos(50))) r 7 . 

Hence, when the conditions on the parameters hold we have 

M{r,6,k,w) < r 4 * aj6 (r) < 

for all r > and we can apply Theorem El In this case since m + = 2 we know 
that system ( 13. 3D has at most two limit cycles, and whenever they exist, one is 
inside the disc T> and is hyperbolic and stable and the other one is outside the 
disc, and is hyperbolic and unstable. Since system (13. 3p has several critical points 
in the disc T>, (0, 0) and (0, ±1) among them, the unstable limit cycle, when it 
exists, surrounds these points. 

An example of parameters for which ty a> b is negative is a = 1/20 and b = 1/15. 
Using the same tools that in the previous case we can prove that the system has 
exactly five critical points and all them are inside T>. Hence by studying the flow 
on the boundary of T> and on {x 2 + y 2 = R 2 }, for R big enough, we prove the 
existence of an unstable limit cycle surrounding T>. 

In short, for these values of the parameters we have proved that system (13. 3p 
has at most two limit cycles and the existence of a limit cycle, which surrounds the 
five real critical points of the system, which are in T>. Our numerical explorations 
indicate that this limit cycle is unique. 

We want to stress that there are very few results in the literature giving upper 
bounds for the number of limit cycles surrounding several critical points. Notice 
that the critical points different of the origin, surrounded by the limit cycle, come 
from the continua of critical points {x 2 + y 2 = 1} that system ( 13. 3ft possesses when 
a = b = 0. The limit cycle is born in {x 2 + y 2 = 2}. 



3.4. Example 4. Consider the system 

x = x(l- (x 2 + y 2 ))(2 - (x 2 + y 2 ))(3 - (x 2 + y 2 )) - y(2 - (x 2 + y 2 )) + ax 2 y\ 
y = x(2- (x 2 + y 2 )) + y(l - (x 2 + y 2 ))(2 - (x 2 + y 2 ))(3 - (x 2 + y 2 )). 

Let us prove that when a is such that 

* a ( r ) : = -98 + ( 192 + — |a| ] r 2 + (-144 + 6|a|) r 4 + (48 + -|a| ) r 6 - 6r 8 < 



for all r > 0, it has at most 3 limit cycles. Once more we use Theorem [A] In 
this occasion we take k = 1, p(s) = (1 — s)(2 — s)(3 — s) and w(r) = r 2 p'{r 2 ) = 
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r 2 (-ll + 12r 2 - 3r 4 ) = -3r 2 (r 2 - 2 - ^§)(r 2 - 2 + Then 

M(r, 0, fc, w) = -98r 4 + — (3072 + 55a sin(20) - 44a sin(40) + 11a sin(60)) r 6 

16 

+ - (-288 - 3a sin(20) + 6a sin(40) - 3a sin(60)) r 8 

+ — (256 - a sin(20) - 4a sin(40) + 3a sin(60)) r 10 - 6r 12 . 
16 

Since m + = 3 and M(r,9,k,w) < r 4 \I/ a (r) < for all r > 0, the existence of at 
most 3 limit cycles follows. For instance the above hypothesis holds for a = 1/34. 
For this value of the parameter we prove, by using the same tools that in the 
previous examples, that the system has several critical points and that, apart from 
the origin, all of them are contained in the ring C = {2 — ^ < x 2 + y 2 < 2 + 
Finally, again similarly that in the previous cases, we prove that there is exactly 
one hyperbolic and stable limit cycle inside the disc {x 2 +y 2 < 2 — ^} and another 

hyperbolic and also stable limit cycle outside the disc {x 2 + y 2 > 2 + ^}. 

The novelty of this example is the existence of a non-trivial system for which 
the maximum number of limit cycles is known (it is 3). Moreover it has at least 
two hyperbolic limit cycles, one of them surrounding only the origin and the other 
one surrounding the first limit cycle and having several critical points between 
them. 

We want to comment that our numerical exploration seems to indicate that 
there is no limit cycle contained in C, and so that the maximum number of limit 
cycles of the system for this value of the parameter is two. 

3.5. Example 5. This last example is interesting from a theoretical point of view. 
Consider the system 

x = xu{x 2 + y 2 ) — yv{x 2 + y 2 ) + eP(x, y), 

y= xv (x 2 + y 2 ) + yu(x 2 + y 2 ) + eQ(x, y), (3.4) 

where u and v are given real polynomials of degree j and assume u is such that all 
their roots are real and simple. Let m + — 1 denote the number of positive real roots 
of u'. Then, for any couple of polynomials P(x,y) and Q(x,y) whose monomials 
have degrees between 2 and 2j + 1, both included, there exists Eo = £q(P, Q) > 
such that if |e| < Eo then system (13.41) is under the hypotheses of Theorem lAl 
Moreover, under these conditions, it has at most m + limit cycles and all the 
existing limit cycles are hyperbolic. 

When e = 0, the above assertions follow from Lemma 12.41 and Theorem [XI 
Consider p(s) = u(s). By using Proposition 12.51 we can take w(r) = r 2 u'(r 2 ) and 
fix a value k > for which 

Pk, w {r) := 2r 4 (p(r 2 )p" (r 2 ) - k(p\r 2 )) 2 ) + 2(1 - k)r 2 p(r 2 )p \r 2 ) < 

for all r G (0,oo). The function p(s,e) for system ( 13. 4p defined in Theorem lAl 
writes as p(s, e) = u(s) + £p(s), for some new polynomial p, also of degree j. It is 
clear that for |e| small enough, the polynomial u(s) + £p(s) has also all its roots 
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real and simple and the number of positive roots of u'(s) + ep'(s) is m + — 1. By 
taking the same value of k and w(r,e) = r 2 (u'(r 2 ) + ep'(r 2 )) we get that the 
function M in the theorem writes as 

dw(r, e) 



M(r, 6, k, w(r, e),e): = R(r, 6, e)- 



dr 



=: Mi{9, k, sY = p k)W {r) + e £ M t (0, k, e)r\ 

i=2 i=2 

where the functions Mj are smooth and 27r-periodic in 9. Since Pk, w { r ) is a negative 
polynomial in (0, oo) of the form Pk,w(f) — &2^ 2 + &4^ 4 + " • •+&4 ? -r 4 - ? , for some b{ real 
numbers and, moreover 62 < and fry < 0, it is clear that for |e| small enough the 
function $fc(r, e) given in Theorem lAl is negative as we wanted to prove. Moreover 
it follows that the maximum number of limit cycles of system (13.4D is m + , and 
whenever they exist they are hyperbolic. 

Notice that it is always true that when a planar system has m + hyperbolic limit 
cycles any small perturbation also has at least the same number of hyperbolic 
limit cycles. The point of the above example is not only to prove the existence of 
at most this number of limit cycles, but to prove that there are planar polynomial 
systems under the hypotheses of Theorem |A] for which all planar systems of the 
same degree near them are also under the hypotheses of the theorem. Moreover, 
as we can see in the study of the previous examples, our approach allows to get 
explicit bounds for the size of the perturbation under which our theorem can be 
applied. 
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